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Abstract 

Semiclassical methods form a bridge between classical systems and 
their quantum counterparts. An interesting phenomenon discovered 
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q h i in this connection is the scar effect, whereby energy eigenstates dis- 

play enhancement structures resembling the path of unstable periodic 
orbits. 

This paper deals with collision states in charged three-body prob- 
ed lems, in periodic media, which are scarred by unstable classical orbits. 

The scar effect has a potential for practical applications because orbits 
corresponding to zero measure classical configurations may be reached 
j_j ■ and stabilized by resonant excitation. It may be used, for example, to 



induce reactions that are favoured by unstable configurations. 



1 Introduction 

Classical and quantum mechanics have qualitatively different features. Nev- 
ertheless, semiclassical methods H provide, in some restricted domains, a 
bridge between the quantum and classical realms. The Van Vleck-Gutzwiller[|2| 
H propagator and the trace formula^] are landmarks in this connection. 

The trace formula relates the fluctuating part of the quantum density 
of states to an oscillatory sum of exponentials, each term corresponding to 
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a classical periodic orbit or to one of its multiple retracings. The actual 
convergence of the periodic orbit sum to the quantum spectrum is still an 
open quest ion 0. Nevertheless many interesting results have been obtained 
by extracting quantum eigenvalue information from the classical periodic 
orbits ||. 

In the trace formula one must sum over all periodic orbits to find a single 
eigenvalue and in Van Vleck-Gutzwiller propagator G(q, q ,t), the momentum 
uncertainty of the q states implies that trajectories of all energies must be 
taken into account. This led some authors to suggest that the reliability of 
the semiclassical results would be improved if one propagates smooth square- 
integrable wave functions with finite-energy uncertainty. The underlying idea 
is that the phase-space localized wave functions would filter out the relevant 
information from the Green's function. 

The most interesting wave function information obtained from semiclassi- 
cal considerations is probably the scar effect^ . For complex systems having 
phase-space domains with sensitive dependence to initial conditions, periodic 
orbits are unstable and even when dense they are nevertheless a zero measure 
set in the smooth ergodic measure. Naively, one would then expect a typical 
wave function to receive contributions from many orbits and its intensity to 
reflect the statistical average of them all. If however the intensity of the 
wave function happens to be either concentrated (or abnormally weak) near 
a classical periodic orbit one says that the quantum state is scarred by the 
periodic orbit. The occurrence of this behavior is easily understood from a 
wavepacket propagation argument || (see below). It has been observed in 
numerical calculations of quantum eigenstates of several classically chaotic 
systems [|| |§ JT(| JTTJ] and was related to unstable periodic orbits through 
semiclassical path integrals [12| |13| fl"4| . 



The scar effect may have far-reaching implications for the practical appli- 
cations of quantum systems. Namely, the unstable classical orbits in chaotic 
systems , even if dense in phase space, are in practice never observed be- 
cause all typical motions are aperiodic and uniformly reproduce the sup- 
port of some diffuse invariant measure. By contrast, in quantum mechanics, 
whenever an unstable periodic orbit scars a quantum energy eigenstate, the 
system may easily be made to behave like the unstable orbit by resonant 
excitation to the corresponding energy level. In this sense scars are a gift of 
Nature, because they allow the exploration of dynamical configurations that 
in classical mechanics are washed away by ergodicity. An interesting scar 
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effect has recently been observed ||15|| on a semiconductor quantum-well tun- 
neling experiment where, by localizing the probability density, the scarring 
of the quantum well states increases the overlap with the emitter and leads 
to enhanced tunnelling for some voltage values. 

In this paper I will be concerned with the relevance of the scar effect in 
accessing collision or near-collision configurations of charged particles in pe- 
riodic media. Of particular interest are scars associated to separatrix orbits 
of which a special type, here called the saddle point scar, is an example. A 
simple wave-packet propagation argument || qualitatively explains the scar 
effect. It also allows the derivation of simple conditions which will be impor- 
tant in the interpretation of the scar effects treated in this paper. Consider 
the overlap integral 

C{t) = {9(t,x)\9(0,x)) (1) 

for a propagating wave packet which at time zero has a Gaussian shape 
and initial conditions (p , x ) corresponding to an unstable periodic orbit. 
Expanding *$f(0,x) in energy eigenstates 

t^j^^W (2) 

n 

one sees that the Fourier transform S(E) of the overlap C(t) is the spectral 
density weighted by the probabilities \c n \ . 

S(E) = Y / \c n \ 2 5(E-E n ) (3) 

n 

On the other hand if the period r of the classical periodic orbit and the 
largest positive Lyapunov exponent A are such that e -rA / 2 is not too small, 
the overlap C(t) will display peaks at times nr. As the wave packet spreads, 
the amplitude of the peaks decreases after each orbit traversal at the rate 
e -rA/2 r-Q^g p our j er transform of C(t) will therefore have peaks of width A 
with spacing uj = — . From Eq.(pJ) one concludes that only the eigenstates 
that lie under the peaks contribute to the expansion of the wave packet. Since 
the wave packet has an enhanced intensity along the region of the period 
orbit, this is expected to carry over to the contributing energy eigenstates. 
This is the scar effect. The stronger the overlap resurgences are, the stronger 
the effect is expected to be. Therefore the intensity of the effect varies like 
1/rA . 
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This simple derivation || of the scar effect is however flawed if the product 
Xd(E) (where d(E) is the mean level density) is very large. Then the num- 
ber of contributing eigenstates is very large and no individual eigenstate is 
required to show a significant intensity enhancement near the periodic orbit. 
Also the argument assumes the low period unstable orbits to be isolated. 
With many nearby orbits having different periods the argument also breaks 
down. Conversely, if it happens that in the same configuration space region 
many different periodic orbits coexist, having the same period, the effect will 
be enhanced. This is the situation for periodic motions in the neighborhood 
of an unstable critical point (a saddle) of a potential function V(x) . Generi- 
cally, in the neighborhood of a critical point, there are coordinates where the 
function may be written as 



On the neighborhood of the point Xi = there are harmonic periodic motions 
along the stable directions (positive cr's) of the critical point. As long as 
anharmonic corrections are unimportant all the orbits have the same period 
independently of their amplitude. The instability parameter of these orbits is 
the smallest negative a (A = — |<7 min |). For each positive <jj the scar intensity 
factor will be -K = j-,/^ , 1 , , where a is the effective mass. The motion 

tA 2tt\ ^ |o- min | ' r" 

in the neighborhood of these unstable periodic orbits being harmonic, the 
energy of the strongest scar is estimated to be 



where the last term corresponds to the ground state energy of the harmonic 
motion on the stable direction and E\ oc is the energy penalty corresponding 
to the localization energy in the transversal directions. In favorable condi- 
tions, that is, if the quadratic approximation for the potential holds over 
a sufficiently large range, higher energy scars might also be observed cor- 
responding to the excited states of this oscillator. They would have the 
approximate energy 




(4) 
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In the following two sections we study a one-dimensional and a three- 
dimensional problem where unstable classical motions indeed leave their trace 
in the some of the quantum states. Of particular interest here are the unstable 
classical motions which correspond to collision or almost-collision configura- 
tions. Some of the potential practical applications of this effect are discussed 
in the conclusions. 

2 Scarred collision states in a one-dimensional 
periodic problem 

Let a quantum system be described by the Hamiltonian 

Wi = -£ - h ~ % + 1/(1x1 - X2l) - n]Xi - y]) - v(lX2 - vl) (7) 

with periodic boundary conditions on a lattice of lattice size L. The Hamil- 
tonian represents a many-body system with two heavy and one light particle 
(7 of order 1CT 4 , for example) in each lattice cell and periodic boundary 
conditions or, alternatively, a three-particle system living on the circle. The 
heavy particles repel each other and attract the light one. For definiteness 
I have considered the heavy particles to be positively charged and the light 
particle to be the negative one. The interaction potential is 

nx) = ^(l + cos(27r|)) (8) 

The factor j- is included in the coupling constant ^ as a convenient factor- 
ization in case one wants to insure scaling properties similar to the three- 
dimensional Coulomb problem. Periodic boundary conditions are imposed 
on the quantum problem by choosing a basis of box-normalized momen- 
tum eigenstates with periodic boundary conditions -j= exp(i2nkx) , k = j 
, n = 0, ±1, ±2, • • •. Denote by\riin 2 p) a state with momenta ^ , ^ and f 
respectively. The matrix elements of Hi in this basis are 

(nVap'lifilmnap) = 5^8^ 5 p/p (2ir) 2 ± (n\ + n\ + ^) ^ 

+ j : 5pp5(n 1 - n[ + n 2 - n^/ifai _ n 'i ~ n 2 + n 2 ) 

-L 6 n' 2 n 2 6 ( n l- n 'l+P-P)fl( n l- n 'l-P + P) 

-L 6 n' ini 6 (P -P +n 2 - n 2 )t\{j> -p -n 2 + n 2 ) 

(9) 
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where fi is the function 

h{a) = J (25a,o + ^,2 + ^,-2) (10) 
The center of mass motion is separated by changing coordinates to 

R = 2^( X l+ X 2 + IV) 

r = X\ — %2 (11) 

v = y-^ 

Classically the dynamics in the center of mass of the three particles is ruled 
by the Hamilton equations 



r = Ap T 



Pr 

5 = 



9V'(r) + lgV'(l- V )-yv'(-^- V ) 



That is, the center of mass dynamics is equivalent to the motion of a two 
dimensional particle in the potential 

U(r, ij) = g |v(r) - v(- - rj) - V(- T - - ry)} (13) 

with an effective Hamiltonian 

H CM = 2p 2 r + ^P^pl + U(r, 77) (14) 



The potential U(r, rj) is displayed in Fig.l with r and rj in units of L. It has 
two stable minima at (r = -|, 77 = 0) (r = — ^, r\ — 0) and a saddle point at 
(r = 77 = 0). The minima correspond to configurations with the two positive 
particles well separated and the negative one midway between the other two. 
The saddle point is a collision state of the two heavy positive particles which 
however has relatively low energy because the repulsive energy of the positive 
particles is compensated by the attractive interaction with the negative one. 
In addition to the unstable fixed point at r = rj = there is a whole collection 
of unstable periodic orbits along the r = axis. According to (H), the energy 
of the corresponding saddle point scar is estimated to be 




E s ~ + E loc + -J » \' u (15) 
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In case one wants to use scar effects to induce collisions the important quan- 
tity to know is the difference between E s and the ground state energy. The 
structure of the ground state will depend on the order of magnitude of the 
physical parameters. For definiteness I assume the ratio of negative to the 
positive particles mass to be very small (O(10~ 4 )) and also of order L -1 in 
the natural units used to write Eq.(|7|). For the numerical results shown below 
the values used are 7 = 2.7 x 10~ 4 , L = 13039 and g = 6. 

From (§) one sees that to excite the first kinetic mode of the light particle 
one needs an energy as compared to kinetic energies of order ^jt- needed 
to localize the heavy particles. Therefore one expects the ground state to 
have a completely delocalized light particle and a heavy particle relative 
coordinate localized around ^ . For this configuration both the transversal 
localization energy and the potential energy are similar to those for the scar 
in (0). Therefore 



A numerical diagonalization of the scaled LH\ Hamiltonian of Eq. (0) was 
performed on a basis of 729 states restricted to zero total momentum (center 
of mass configurations). The structure of the energy spectrum is shown in 
Fig. 2. The band structure corresponds to the energies needed to excite the 
successively higher kinetic modes of the light particle. The first state in each 
band lies approximately at the energy 



n = 0, 1, 2, • • • . The ground state energy is E Q = —5.91. The amplitude of 
the ground state wave function is shown in Fig. 3. The amplitude of the state 
at the top of the negative energy band (E s = —0.905) is shown in Figs.4a,b. 
This is the state that corresponds to the saddle point scar described above. 
The difference E s — E is 5.01 to be compared with the value 5.8 obtained 
from the analytical estimate (|ToT). The state, at the top of the first band, is, 
in this band, the one with the largest heavy particle overlap, defined as 




(16) 
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For other bands it also happens that the state with the largest overlap is at 
the top of the band. I do not know whether there is a general mechanism 
leading to this fact or whether it is a particular feature of this model. There 
are in each band other states with large overlaps which however, as a rule, 
are not so concentrated along the r = line (see for example in Fig. 5 the 
contour plot for the state at E — —1.298 ). 

In this system, the harmonic approximation used in (||) to predict the 
first scar state energy, cannot be reliably carried over to higher harmonic 
excitations, because higher energy states are not so localized around (r = 
f] = 0) and the quadratic approximation for the potential is no longer valid. 
For example for the state at the top of the second band (Fig. 6) we have 
E — Eq = 18.67 whereas ([]) with n = 1 would predict 17.4. For the state at 
the top of the third band (Fig. 7), E — E = 51.22 whereas @ with n = 2 
yields 29.0. 



3 A three-dimensional periodic Coulomb prob- 
lem 

The three-dimensional problem considered here is, as before, a periodic bound- 
ary condition many-body problem with two heavy and one light particle in 
each lattice cell of volume L 3 or, alternatively, a three-particle system living 
on the 3-torus. The interaction is the Coulomb potential 

V(\x 1 -x 2 \) = -. (17) 

\X\ — x 2 \ 

In problems with periodic boundary conditions and long-range interactions 
it would make sense to consider that, acting on the particles in one cell, are 
the forces of all the particles in the other cells or at least of those on the 
neighboring cells. Alternatively one might change the potential to make it 
L— periodic. 

Here periodicity is explicitly introduced by the choice of a periodic basis 
of momentum eigenstates as in the one- dimensional problem. However for 
the computation of the matrix elements of the Hamiltonian only the forces 
between the three particles will be considered. This slightly underestimates 
the interaction part of the Hamiltonian but it does not change the qualitative 
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nature of the problem especially in what concerns the small distance collision 
effects. 

For this problem I will first find the low-lying quantum spectrum and 
then see whether the strong overlap states may or may not be interpreted 
as scar states. In a basis of momentum eigenstates \nln 2 rf) with \n) = 
(L)~2 exp(i^ 7? .it ), the matrix elements of the Hamiltonian 

H 2 = -Ai - A 2 - -A 3 + \ Xl - x 2 \- 1 - \x x - y^ 1 - \x 2 - y^ 1 (18) 

7 

are 



n 1 n 2 p \tL 2 \nin 2 p 



n x n\ n 2 n 2 p p 
Pi? 

n 2 r? 2 _^ 



nif + \n 2 \ 


2 + M! 


) 


~n 2 )f2{\n[ 


- 


n 2 + n 2 |) 


P )/2(K 


- w i - 


P + P 1) 




- p - 


n 2 +n 2 |; 



(19) 



with 



/ 2 (H) = — (l-oos (ttpH)) (20) 



and p = 2 ( 3 . For the computation of the spectrum of H 2 one considers 
also states symmetric and antisymmetric under interchange of the two heavy 
particles. 

~ft)± = - (|nin2 ~t) ± K^2^)) (21) 
c 

with c = for nt 7^ n| and c = 2 for nt = n|. For the numerical diago- 
nalization of the Hamiltonian H 2 a basis of 3176523 states was considered. 
Using momentum conservation and conservation of permutation symmetry 
type the matrix is however reduced into invariant blocks to simplify the 
computation. Fig. 8 shows the lower part of the spectrum for center of mass 
(zero total momentum) states with symmetric states denoted by crosses and 
antisymmetric ones by dots. In Figs. 9a and 9b one compares the lowest ly- 
ing antisymmetric (^oa) an d symmetric (^os) states. The quantity that is 
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plotted is the integrated probability 

l*(M,M)r = / d£l r dn n \*{r,ri)\ 2 (22) 

The lowest lying symmetric eigenstate displays a strong overlap of the heavy 
particles. The nature of this state is better understood from the two-coordinate 
projections 

|^( ri ,^)| 2 = J d 2 rd 2 r]\m (r,r/)| 2 (23) 

the integration being carried over all coordinates other than r, and rjj. Figs. 10a 
and 10b show |\& (r 1; r/i)| 2 and \^ {ri,r] 2 )\ 2 ■ All other like-coordinate and 
unlike-coordinate projections are identical to those shown in the figures. 

As in the one-dimensional case studied before, there is a simple relation 
between this state and classical orbits. The classical equations for center of 
mass motion are 




(24) 

For collisions or near-collisions to take place at low energies the light parti- 
cle must be near the two heavy particles part of the time in order for the 
attractive energy to compensate the strong repulsion. This implies some 
localization of the light particle in one dimension at least. Localization 
in more dimensions however is very costly in kinetic energy. Therefore it 
is likely to find, associated to the low lying symmetric state, orbits corre- 
sponding to planar motion but not to collinear motion. From ( p4|) it fol- 
lows that there are orbits in the plane (7*1,772), that is, if at t = one has 
7*2 = Pr2 = r 3 = Pr3 = Vi = Pr)\ = V3 = Prfi, = then the same holds true for 
all t. The phase-portrait in the plane (7*1,772) is symmetric about the 7*1 = 
line. In Fig. 11a and lib I have plotted two typical orbits. Periodic boundary 
conditions are imposed at | and — |. In-between the orbits that move to the 
right and to the left there is the separatrix at T\ — that passes through the 
singular point of the potential. It is precisely along this separatrix that the 
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state \l/os is strongly scarred. In this case we have not found a saddle scar 
in the sense defined in the introduction because the separatrix goes through 
a singular point of the potential. It is possible that by a change of time, 
and using the well-known regularization|l6[ of the Kepler problem one may 
still be able to formally interpret the state as a saddle scar. This is not 
however very important. What is important to notice is that once again we 
have found that the low-energy collisional state is associated to an unstable 
feature of the classical phase-portrait. 



4 Conclusions 

Through the scar effect, orbits that correspond to zero measure classical con- 
figurations may be reached and stabilized in quantum mechanics by resonant 
excitation with the appropriate energy. Their location in the energy spec- 
trum may, in favorable cases, be found either from theoretical considerations 
(from symmetry, from being at the top of bands, etc.) or from absorption 
experiments. 

This characteristically quantum phenomenon, may be practically used to 
induce reactions which are favored by unstable or difficult to reach configura- 
tions. An obvious potential application is to fusion reactions. Most practical 
nuclear fusion mechanisms proposed so far have a two-mechanism nature. 
For example muon catalyzed fusion relies on the fact that the muon is 200 
times heavier than the electron to bring the bound nuclei together, but then 
it is the tunneling effect than will eventually allow the nuclei to fuse. In 
toroidal plasma confinement it is the magnetic field configuration that keeps 
the ions together, but then it is radio frequency or ion injection heating that 
supplies them with enough kinetic energy to overcome the Coulomb barrier. 
What I am proposing here is that if enough deuterons, for example, are con- 
fined in a periodic medium (a crystal lattice, for example) then, resonant 
excitation and the scar effect may be used to excite collisional states and 
induce fusion reactions. I strongly emphasize that like in the known fusion 
methods we should separate the two problems of confinement and collision. 
The crystal lattice only serves as a confinement device, an additional colli- 
sion mechanism being needed, which the scar effect, discussed in this paper, 
may provide. This is contrary to the hopes of the cold fusion saga where 
spontaneous fusion was expected just from confining deuterons in a lattice. 
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Actually a simple calculation shows that in the problem discussed in Section 
3 the separation between the ground state and the first scar state is such 
that thermal excitation is highly improbable. However resonant excitation 
by electromagnetic radiation seems possible. 

For other applications of the scar effect and in particular of its saddle 
point enhancement one might think of catalyzing chemical reactions on lattice 
substrates. 

In conclusion: unstable configurations corresponding to unstable periodic 
orbits are in classical mechanics of little use unless very sophisticated control 
methods are usedJCT]. This is because, for each energy level, the stable 
invariant measures are smoothly distributed all over the energy surface. In 
this sense quantum mechanics is more parsimonious because for each energy 
level it displays just a fraction of the blurred picture of smooth classical 
dynamics. In particular, by isolating through the scar effect improbable 
classical configurations, quantum mechanics is, for practical applications, a 
cure to the classical mechanics curse of ergodicity on the energy surface. 

5 Figure captions 

Fig.l - Effective center of mass potential for the one-dimensional problem 
Fig. 2 - Energy spectrum in the one-dimensional problem 
Fig. 3 - Ground state wave function 

Figs.4a,b - Wave function of the state at the top of the negative energy 
band 

Fig. 5 - Contour plot for the state at energy E = —1.298 
Fig. 6 - Wave function of the state at the top of the second energy band 
Fig. 7 - Wave function of the state at the top of the third energy band 
Fig. 8 - Energy spectrum in the three-dimensional problem 
Figs.9a,b - |^(|r|, |?7|)| 2 for the lowest lying antisymmetric and symmetric 
states 

Fig. 10a - |\l/(ri, ?7i)| 2 for the lowest lying symmetric state 
Fig. 10b - |\I/(ri, 772) | 2 for the lowest lying symmetric state 
Figs.lla,b - Two classical center of mass orbits in the plane (ri, 7/2) 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/quant-ph/9705040vl 
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